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Accurate Monte Carlo data from a set of isotherms near the critical point are analyzed using two 
RG based complementary representations, given respectively in terms of h=h/\t\ l3S and f=t/h}^ s . 
Scaling plots for data on simple cubic Ising lattices are compared with plots of ML 13 " vs. \t\L}l v 
for increasing L values and with high quality experimental data on C'rBr^. Finite size effects and 
the equation of state are discussed. 
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Monte Carlo methods using Wolff algorithms 
have been extensively used to describe M(T) for H=0 in 
the whole range of temperatures both below and above 
the critical temperature (specially at T=T C ) but few, 
if any, Monte Carlo simulations of magnetic isotherms 
M(H) at T=T C have been reported in the literature. 
It is clear, however, that with recent improvements in 
computing facilities (larger memory, greater speed, bet- 
ter availability) accurate, well thermalized, closely spaced 
data can provide very substantial contribution to describ- 
ing the phase transition and better understanding of fi- 
nite size effects. 

Simulations of this type are reported in this work, per- 
formed using Metropolis algorithms 3], which are spe- 
cially convenient to describe the system evolution at con- 
stant temperature and for small field increments^, 0. 
Using relatively large lattices with 70x70x70 spins, an 
accurate characterization of the scaling behavior, and, 
therefore, the equation of state in the vicinity of the 
critical point (H=0, T C =4.5II523785) 6] for simple cu- 
bic Ising lattices can be obtained. Details of the Monte 
Carlo calculations for isotherms taken at T near T c are 
given in reference . Periodic boundary conditions were 
used and 140,000 Monte Carlo steps were taken at each 
field/temperature to ensure equilibrium. 

For the scaling representation of the raw M(H) data at 
each T we did use two complementary ways: (a) the usual 
way 0, involving fh=M(t, h)/\t\P and h=h/\t\P s , where 
t = (T—T c ) /T c involves the temperature gap to T c , h (re- 
duced field) is proportional to H, (3 is the spontaneous 
magnetization critical exponent (3 = dlagM / dlog\t\ at 
T^T C , and S is the critical isotherm exponent 5 = 
91ogM/<91og/i at T-^T C ; and (b) a complementary way 
involving p,=M(h, tj/h 1 / 5 and f^t/h 1 / 135 . 

Near the fixed point corresponding to the critical point, 
the singular part of the reduced free energy per spin may 



be written in scaling form as 

f s ( 9l , g 2 , ga, . . .) ~ b d f s (b^ 9l ,b^g 2 , b^g 3 , . . . ) (I) 

where b is an arbitrary scale factor and the yi relate to 
the usual critical exponents. To get (a) Equation Q is 
differentiated with respect to the field in the usual way 
to obtain 

M(t, h,g 3 , . . . ) ~ b-^MQPH, V»h, 6 y3 33, • • ■ ) (2) 

Taking 6 yi |£| = I and setting the irrelevant variables 
equal to zero 

M(t, h) ~| t M(±I, h | t |-»M) (3) 

and using the well known scaling relationships giv- 
ing yi and y 2 in terms of j3 and 5 one gets 

M(t, h) ~| t f M(±l, h | t I" 133 ) (4) 

which implies scaling using the ordinary scaling vari- 
ables 

m = M(t,h)/ \ tf,h = h/\t f s (5) 

To get (b), on the other hand, we can argue in an anal- 
ogous way. Taking b V2 h=l in Equation (J2J and setting 
the irrelevant variables equal to zero, 

M(t, h) ~ h( d - V2)/v2 M(t ■ h~ Vl/v2 , 1) (6) 

and using y\ and y 2 in terms of f3 and S, one finally 
gets 

M(t,h) ~ /i 1 / 5 M(*- 1//3,5 , 1) (7) 
which implies as alternative scaling variables 
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tL = M{t.K)/h 1 l & ,T = t/h 1 l f3S (8) 
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FIG. 1: 

lattice of 70 
fi=M/h 1/s 



Scaling plots of Monte Carlo data for a s.c. 
3 spins, (a) 
t/h 1 ^ 6 . 



Ising 



m=M/\tf vs h=h/\tf s , and (b) 
vs T=t//i"" J ". The data include 30 isotherms in 
the intervals 4 < T < 5 and < h < 0.02, the same symbol 
has been used for all of them. The critical temperature I™ 
was taken as 7^=4.511523785 and the critical exponents p| 
as ,0=5/16 and 8=5. 



Figure 1(a) gives our Monte Carlo data using m and 
h as scaling variables as given by Equation JSJ with 
T c =4.511523785 and the exponents /3=5/16=0.3125 and 
5=5. The usual behavior is observed. We may note that 
the data scale extremely well and that only very minor 
deviations at lower h, attributable to finite size effects are 
perceptible. Figure 1(b) gives the same data using the 
alternative scaling representation. It may be noted im- 
mediately that the plot in Figure 1(b) resembles closely 
typical scaling plots of ML 13 /" vs \t\L x l v for lattices 
with linear dimension L, in our case, for Ising systems 
of L 3 spins, suggesting [8| formal relationships between 
H and L, M and L, and \t\ and L, which, through the 
scaling plot branches defining the critical isotherm, the 
spontaneous magnetization (coexistence curve) below T c , 
and the low field susceptibility above, imply respectively 




FIG. 2: Scaling plot of Monte Carlo data ML W " vs \t\L x l v 
for s.c. Ising lattices with linear size L=30, 60, 90, 115. Note 
that finite size effects for L=30 show up closer to \t\— »0. 



TABLE I: Critical exponents 





P 


8- 1 


7 


Ising 3D [4] 


5/16=0.3125 


1/5=0.2 


5/4=1.25 


Heisenberg 3D [7] 


0.340 


0.208 


1.39 


CrBr 3 [9] 


0.368 


0.233 


1.215 



H ~ L'W", M ~ L~ l v and \ t |~ L~ 



l/u 



(9) 



Hence 

ML / V ^const (i^0) — ► critical isotherm 
ML /" ~constx(|£|L 1 / l ') /3 (t < 0) — > spontaneous 
magnetization 

ML /" ~constx{\t\L 1 /»y 3(5 - 1 '>/ 2 (t > 0) — ► suscepti- 
bility 

Figure H shows data of ML /" vs \t\L x l v for H=0 
and L=30, 60, 90, 115 which mimic the behavior shown 
in Figure 1(b) implying that simulations using periodic 
boundary conditions of phase transitions with finite size 
show the effects of an effective straining contribution to 
the magnetic field Hf s ~L~~ 0S / U , i.e. H e ff=H+Hf s such 
that for L— >oo, H e ff=H. 

Figure 3 gives scaling plots of the high quality data of 
30 isotherms at the vicinity of the Curie temperature per- 
taining to the insulating ferromagnet CrBr^ measured by 
Ho and Litster 0, which are a classical example of ex- 
perimental scaling data. We show plots of m vs h (Figure 
3a) and p, vs f (Figure 3b) with T = T C =32.844K for two 
sets of critical exponents, Ising 3d (fractional values) and 
Ho & Litster (experimental values) summarized in Table 
HJ Both sets of critical exponents produce good scaling 
plots in (a) as well as in (b), but the critical exponents 
of Ho and Litster produce somewhat better scaling plots. 
The accuracy of the magnetization measurements for the 
set of isotherms was comparable to that of nuclear mag- 
netic resonance data and it was sufficiently precise to es- 
tablish the form of the scaling function. Our Monte Carlo 
data, shown in Figures 1(a) and 1(b) are of comparable 
quality to establish the 3d Ising scaling function. 
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FIG. 3: Scaling plots of experimental data for CrBr^ (a) 



M/\tf vs h=h/\tf, and (b) p,=M/h 1/d vs f=t//i 



The data are made up of 30 isotherms in the interval T C -0.9K 
< T < T C +6.7K. The critical temperature was T c =32.844 K 
and the critical exponents used were (3=5/16, 8=5 (Ising 3D) 
and /3=0.368, 5=4.28, as given in Reference [Sj. 
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Finally we address the question of the form equation 
of state for 3D Ising lattices in the light of the infor- 
mation provided by the set of isotherms in the vicinity 
of the critical temperature obtained by the Monte Carlo 
method in our 70 3 s.c. lattice. Figure 4 gives the plot of 
M(h,t), Equation (JSJ), rewritten as 



f(m) = A{l±Bm z ) 



(10) 



Here A can be reduced to unity just by choosing prop- 
erly the units for the field H. B is a more meaningful co- 
efficient which, in the particular case of a phase transition 
describable by means of the mean field approximation 
(such as the phase transition in a uniaxial ferroelectric) 
is equal to 1/5=1/3. And z, as pointed out in Reference 
is given by z=05/v=2.5 for T <T C and z=P5=\M2 
for T > T c . Figure 4 shows the excellent fit obtained by 
means of Equation JTUJ with (f?/A)=0.102. The equa- 
tion of state put in the form given by Equation IjlUfl is 
specially good to show directly the most relevant infor- 
mation: (a) the critical isotherm for T < T c and T > T c , 
(b) the spontaneous magnetization curve (T < T c ) as a 
vertical line, and (c) the zero field susceptibility (T > T c ) 
as a horizontal line. The quality of the fit is comparable 
or better than those obtained with traditional ex pres sions 
of the scaling function @, IH 13 13 III 13 13 113 • 

Work is in progress to obtain Monte Carlo data in 
larger 3D Ising lattices at closer field/temperature in- 
tervals, and to extend the investigation of scaling plots 
in the vicinity of the transition to higher dimensionali- 
ties Ising 4D, Ising 5D, etc, in order to monitor closely 
how the approach to mean field behavior takes place. Of 
course we will be limited to more reduced sizes (smaller 
L's) as the dimensionality increases, but we have excel- 
lent experimental data |18( on a complete set of isotherms 
in uniaxial ferroelectric TGS to produce excellent scaling 
plots with T c =321.470 and mean field critical exponents 
/3=l/2 and 8=3. 



Acknowledgments 

We specially acknowledge helpful comments and soft- 
ware by M.I. Marques. Support from the Spanish 
MECyD through Grant Number BFM2000-0032 is grate- 
fully acknowledged. 



FIG. 4: Optimized fits of the 30 Monte Carlo isotherms to 
the given Ising 3D equation of state. 
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